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Abstract 

The refined enumeration of alternating sign matrices (ASMs) of 
given order fiaving prescribed behavior near one or more of their bound- 
ary edges has been the subject of extensive study, starting with the 
Refined Alternating Sign Matrix Conjecture of Mills- Robbins-Rumsey 
[21j . its proof by Zeilberger [27], and more recent work on doubly- 
refined and triply-refined enumeration by several authors. In this paper 
we extend the previously known results on this problem by deriving ex- 
plicit enumeration formulas for the "top-left-bottom" (triply-refined) 
and "top-left-bottom-right" (quadruply-refined) enumerations. The 
latter case solves the problem of computing the full boundary cor- 
relation function for ASMs. The enumeration formulas are proved by 
deriving new representations, which are of independent interest, for the 
partition function of the square ice model with domain wall boundary 
conditions at the "combinatorial point" rj = 27r/3. 

1 Introduction 

1.1 Alternating sign matrices 

An alternating sign matrix (ASM) of order n is an n x n square ma- 
trix with entries in {0,1,-1} such that in every row and every column, 
the sum of the entries is 1 and the non-zero terms appear with alternating 
signs; see Figure [T] for an example. From this seemingly innocuous defini- 
tion a uniquely fascinating class of objects arises: originally discovered by 
Mills, Robbins and Rumsey in connection with their study of Dodgson's 
condensation method for computing determinants, ASMs have since been 
found to have deep connections to many other topics of interest in combina- 
torics and statistical physics. Some places where ASMs make an unexpected 
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Figure 1: An alternating sign matrix of order 6. 



appearance are: the square ice model (a.k.a. the six-vertex model) 
[ini [251 [27] ; totally symmetric self complementary plane parti- 
tions [I6l[22]; descending plane partitions [21]; domino tilings fll] : 
and the 0(1) loop model in a cyhndrical geometry [9l I24j . 

1.2 Enumeration of alternating sign matrices 

The current paper will be focused on one particular aspect of the study of 
ASMs, namely the problem of enumerating various naturally-occurring sets 
of ASMs of some fixed order n. This is a problem with a venerable history, 
starting with the seminal paper [21] of Mills, Robbins and Rumsey. Having 
defined ASMs and discovered the role that they play in the definition of 
the A-determinant, a natural generalization of matrix determinants. Mills 
et al. considered the problem of finding the total number An of ASMs of 
order n. Based on numerical observations, they conjectured the formula 



Another natural enumeration problem concerned the so-called refined enu- 
meration of ASMs. It is based on the trivial observation (an immediate 
consequence of the definition of ASMs) that the top row of an ASM con- 
tains a single 1 and no —Is. The position of the 1 in the top row is therefore 
an interesting parameter by which one may refine the total enumeration. 
Thus, for 1 < k < n Mills et al. defined 




(1.1) 



^n,k = # of ASMs of order n with 1 in position (1, k) 



and conjectured that 
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The conjectures (1.1) and (1.2) became the subject of intensive study by 
combinatorialists. The former became known as the Alternating Sign Matrix 
Conjecture, and the latter as the Refined Alternating Sign Matrix (RASM) 
Conjecture; both were eventually proved by Zeilberger \26\ 27J (in the case of 



(1.2), building on techniques introduced by Kuperberg in [20j). A readable 



account of these developments can be found in the book by Bressoud [6]. 

1.3 Doubly-refined enumeration and beyond 

By symmetry, the observation mentioned above concerning the behavior of 
the top row of an ASM applies not just to the top row but also to the bottom 
row and to the leftmost and rightmost columns. This lead Mills et al. to 
consider a doubly-refined enumeration involving two parameters, one for the 
position of the 1 in the top row and another for the corresponding position 
in the bottom row. Thus, for 1 < i, j < n we may denote 

^™(i,j) = # of ASMs of order n with 1 in positions (l,i), (ra,j). 

Mills et al. did not conjecture an explicit formula for ^™(i,j). However, 
they discovered that a different family of objects, namely the so-called to- 
tally symmetric self-complementary plane partitions (TSSCPPs), had a two- 
parameter refinement which they conjectured J22J is also given by the same 
family of numbers (A™(i, j))-^^^ This was proved in recent years by 
Fonseca and Zinn- Justin [16] . Also fairly recently, Stroganov [25] derived an 
explicit formula for these "top-bottom" enumeration coefficients, expressing 
them in terms of the singly-refined coefficients An^k- His result states that, 
for 1 < i < j < n, 

i-l 

Al''{i,j)=An,j-i + ^Dn{k,j-i + k), (1.3) 



k=l 



where 



Dn{s,t) = {An-l,t{An,s+l " ^n.s) + ^n-l,s(-4„,t+l - ^„,t)) . 

(To be a bit more precise, Stroganov proved that the numbers An{i,j) sat- 
isfy the recurrence relation An{i + 1, j + 1) — An{i,j) = Dn{i,j) for all 



^ ^ ijj ^n; this immediately implies (1.3) by summation.) Stroganov also 
considered another natural doubly-refined enumeration, the "top-left" enu- 
meration of ASMs with prescribed behavior in the top row and left column. 
For 1 < i, j < n denote as above 

^n^ihj) = of ASMs of order n with 1 in positions (l,i), (j, 1). 



3 



Stroganov proved that the top-left coefficients {A'J^^{i, are related to 
the top-bottom coefficients {A^^{i, via the relations 



= A^^(i,j + i) + A^^(i + i,j)-^n + + (^i >2), 

(1.4) 



It is not difficult to invert these relations and therefore to obtain an explicit 
formula for A'^^{i,j): 

' An-\ ifi = j = l, 

\ii = \<2 

or J = 1 < i, 



p=l 9=1 



(The derivation of this identity from (1.4) involves an elementary manipu 



lation of generating functions and is omitted; however, note that this very 
interesting formula does not seem to have appeared previously in the liter- 
ature.) 

A further development of the theory of refined enumerations of ASMs 
came with the papers [U] , |19] . In [H] , Fischer and Romik introduced a new 
family of doubly-refined enumeration coefficients, which we will call here the 
"top-two" (or "top-top" ) coefficients, since they enumerate ASMs based on 
their behavior in the top two rows. Let 

X #|m = (mij)"j=i e ASM„ I nii^i + m2,i = mij + = l|, 

where we use the self-explanatory notation ASM„ to refer to the set of 
ASMs of order n (the division by j — i + 1 is used to factor out the additional 
condition mi^t = 1 that can be — but isn't — imposed on the matrix for any 
i ^ k < j; see ^4j for more details.) Fischer and Romik studied the problem 
of finding a formula for A]^^{i,j). They managed to derive a system of 
linear equations satisfied for each n by the coefficients {Aj^'^{i,j))ij, which, 
modulo a reasonable conjecture on the invertibility of the system, could 
be used to express them explicitly as ratios of determinants by Cramer's 
rule. They also conjectured a more explicit but very complicated formula 
expressing A'^'^{i, j) in terms of a hypergeometric summation. Finally, in the 
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paper [19], Karklinsky and Romik derived a much simpler explicit formula. 
They proved that, for 1 < i < j < n, 

n-j P /X 
p=0 q=0 

where 

En{s,t) = — {An-l,t{An,s+l — ^ra,s) — ^n-l,s(^n,4+l — ^n,t)) • 

(Note the similarity in the definitions of Dn{s,t) and En{s,t), which differ 
only by a single sign; it is also easy to check that En{s, t) = Dn{s, n — t) by 
using the symmetry An^k = ^n,ra+i-fc-) Another proof of this formula using 
different methods was given shortly afterwards by Fischer |12j . 

Going beyond the doubly refined enumeration, Fischer and Romik started 
developing a theory for the "A;-tuply refined" enumeration of ASMs with pre- 
scribed behavior in the top k rows. For each A; > 1, they defined a family 
of coefficients {A^ {ji, ■ ■ ■ , jk))i<ji<...<jk<n that encode this enumeration 
(the correct way to do it is to think about ASMs as complete mono- 
tone triangles; see I14j), and conjectured a generalization of the system of 
linear equations proved to hold for the top- two doubly-refined coefficients. 
That conjecture was later proved by Fischer [12]. Recently, Fischer also 
derived additional linear equations relating the "top-left-bottom" and "top- 
top-bottom" families of triply- refined enumeration coefficients [13]. How- 
ever, these results left open the problem of deriving explicit, closed- form 
formulas to compute these coefficients. 

1.4 New results: triply- and quadruply-refined enumeration 

Our goal in this paper is to extend the known results on enumeration of 
alternating sign matrices with prescribed behavior at one or more of their 
boundary edges, to fully take into account the joint behavior with respect 
to all edges, thus effectively completing the study of this type of boundary 
enumeration. 

We consider the triply-refined and quadruply-refined enumerations for 
ASMs with prescribed behavior in three (resp. four) of their boundary 
rows/columns. The triply refined enumeration coefficients, which we will 
refer to as "top-left-bottom" coefficients, are given for any 1 < i,j,k < nhy 

Al^^{i,j, k) = #{Af = (mij)^^=i G ASM„ | mi,i = mj,i = = l}. 
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Similarly, the quadruply-refined "top-left-bottom-right" enumeration coeffi- 
cients are defined for 1 < i,j,k,i < n by 

Al^^^{i,j,k,e) = #{(mij)^j=i G ASM„ | mi,i = mj- 1 =m„,fc = m£,„ = l}. 

We remark that in the statistical physics literature, the vector of enumer- 
ation coefficients measuring prescribed behavior along specific boundary 
edges (or in some situations the generating function of this vector) is often 
referred to as the boundary correlation function. Boundary correlation 
functions have been studied in great generality by others in the context of 
integrable systems [H El [121 E3] • 

Our first main result is the following explicit formula relating the top- 
left-bottom coefficients to the usual singly-refined coefficients (^n,fc)fc=i- It 
turns out to be most natural to express the formula as an identity between 
two multivariate generating functions. 

Define (single- variable) generating functions an{t) , (3nit) , jnit) and (5„(t) 

by 

n 

k=l 
n-1 

Mt) = Y,^n-l,kt^ =tan-l{t), 
k=l 
n+2 

ln{t) = ^ - 2(n - + 3)An-l,fc-3 + (5n - 4/c + Q)An-l,k-2 
k=l 

+ (n + 4A; - 6) A„_i,fe_i - 2A;^„_i,fc) 

n+3 

^n{t) = Y^ (A{n + A-k){n + h-k)An-i,k-r, 

k=l 

- 4(n + 4 - /c)(5n + 11 - 4fc)A„_i_A.^4 

+ (240 - 172fc + 32fc^ + 120n - 52/cn + 2W)An-i,k-^ 

- 2(80 - 80A; + 2{)k^ + 42n - 20/cn - f)n^)An-i,k~2 
+ (64 - 84A; + 32/0^ - An - 12kn + n2)A„„i^fe-i 

- ikin - 5 + 4:k)An-i,k + 4:k{k + l)An-i,k+i)t^~\ 

and denote by A(ti, . . . ,tk) = ni<i<j<fc(*« ~ the standard Vandermonde 
product of real numbers ti, . . . ,tk- 
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Theorem 1 (Triply refined boundary correlation function). The generating 
function of ASMs refined according to the positions of the Is in the first row, 
leftmost column and last row satisfies 

(n n—1 n 
i=2 j=2 k=2 

/ {x-lfonix) (y-l)2a„(y) {z - l)^an{z) 

= PnA{x,y,z)-Uetl (x - l)/3„(x) {y-l)Pn{y) (z - l)/3„(z) 

\ Inix) Jn{y) 7niz) 

n 

-{l-y + xy)Y, An-l,^-lX*-'2"-' 

i=2 
n 



k=2 



(1.6) 



where pn is a constant given by 



1 / (n-2)!(3n-2)! y 
^" 8A2(2n-2) V(2n-3)!(2n- l)!y ' ^ ' 

Our next result gives an analogous expansion for the generating function 
of the quadruply refined enumeration coefficients {A]^^^^{i,j,k,£)). Natu- 
rally, the expressions become more complicated, but the result is structurally 
similar. It will be convenient to define the doubly refined generating func- 
tion, 

n 

Al'^{x,y) = ^ Al^{i,j)x^-'y^-'. (1.8) 

Note that A]i^{x,y) can be expressed up to a constant in terms of the 
functions a^, and Pn defined in (1.5) by combining equations (32) and (33) 
of 



Theorem 2 (Full boundary correlation function). The generating function 
of ASMs refined according to the positions of the Is in the first row, leftmost 
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column, last row and rightmost column satisfies 

4 /n— 1 n— 1 n— 1 n— 1 

t=l \ i=2 j=2 k=2 e=2 



CJ„A(xi,X2,X3,X4) ^ 

^ (xi a„{xi) (x2-l)^o„(x2) (x3-1)3q„(x3) {x4-l)^a„{x4) ^ 

X det 



i-2 



{xi-l)^l3„(xi) {X2-1)^I3„{X2) {X3-1)2/3„{X3) (a;4-l)2/3„(x4) 

(xi-l)7„(a;i) (X2-I)7„(x2) (3:^3-l)7n(a;3) (a;4-l)7n(^4) 

\ 5n(a:i) Sn{x2) Sn{x3) Sn(x4) j 

4 2 



^ JJ(1 - xt+fc + xt+k~\xt+k) x1_lx'l^l At 

t=l k=0 

An~2 (1 - 2:3 + X2X3)(1 - Xi+ X4X1) X^'"^ X^'"^ 
An-2 (1 - X2 + XiX2)(l - X4 + X3X4) X^~'^ , 



-1 { xu — ^ 

Xt+1 



(1.9) 



with the convention that xt = xt-4 for 5 < t < 8. The constant an is given 
by 

n 1 / (n-2)!(3n-2)! y 

'^n = ^TTTTTT- 777- ^TTTTT- TVT • (1.10) 



64^3 (2n - 2)2(2n - 3) \{2n - 3)!(2n - 1) 



Note that the identities (1.6) and ( |1.9| ) make it possible to compute 
the coefficients X^^^{i,j,k) and A^^^^{i, j,k,l) in a computationally ef- 
ficient manner using standard polynomial algebra. The Maple package 
Ref inedASM1234 accompanying this paper, which may be downloaded from 
the authors' web pages, provides a demonstration of the practical application 
of these formulas. 

The basis for our enumeration results is a new expression for the partition 
function of the square ice model with the so-called "domain wall" bound- 
ary conditions at a special value of the parameter known as the "crossing 
parameter"; see Theorem [5] below. This formula is of independent interest, 
both for its intrinsic theoretical value and because of its applicability to the 
problem of refined enumeration. In the next section we recall the basic facts 
about the square ice model and develop this new approach to the study of its 
partition function. In section 3 we apply these results to prove Theorems [T] 
andU 
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2 A new expression for the square ice partition 
function 



The square ice model, also known as the six-vertex model, is a well known 
"exactly solvable model" from statistical physics [2|. The connection to 
ASMs appears when the model is considered with specific boundary condi- 
tions. We start with an n x n square lattice with 4n half-edges emanating 
from the boundary vertices. To the ith row, we associate a real-valued pa- 
rameter Xi and similarly, to the jth column, we associate a parameter yj. 
These are collectively called the spectral parameters, with the Xj's re- 
ferred to as the row parameters and the y^-'s termed the column param- 
eters. A state, or configuration, of the model is an assignment of arrows 
to each edge (including the half-edges at the boundary) so that every vertex 
has two incoming and two outgoing edges. There are six types of vertices, 
depicted in Figure [2] below. 




(typel) (type 2) (type 3) (type 4) (type 5) (type 6) 
Figure 2: The six types of vertices in a square ice configuration. 



Type 2 is the rotation by 180° of type 1, type 4 is the rotation by 180° 
of type 3, and type 3 is the vertical reflection of type 1. Types 5 and 6 are 
rather special because in both these cases, the two vertical (as also horizontal 
arrows) both point inwards or both point outwards. The domain wall 
boundary conditions (DWBC) correspond to the horizontal half-edges 
at the boundary pointing inwards and the vertical ones pointing outwards. 
See Figure [3] for an example. 

The relevance of square ice to the study of ASMs comes from the well- 
known fact that configurations of the square ice model with domain wall 
boundary conditions are in bijection with alternating sign matrices. Given 
a square ice configuration satisfying the DWBC, it can be translated to an 
alternating sign matrix by replacing every vertex of types 1 through 4 with 
a zero, vertices of type 5 with — I's and vertices of type 6 with I's. As an 
example, the configuration in Figure [3] maps to the ASM in Figure [T] under 
this bijection. 

Fix six functions Xt(x, y), 1 < t < 6. Given a square ice configuration C 
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Figure 3: An example of a square ice configuration with domain wall bound- 
ary conditions. 



of order n, for each 1 < i, j < n we assign a weight . ■ (xj, yj) to the vertex 
(i, j), where Uj is the type of the vertex at (i, j). We then define the weight 
of the configuration C to be the product of the weights of all the vertices, 

w{C)= n Xt,,{xi,yj). (2.1) 

l<i,j<n 

Letting Cn denote the set of square ice configurations of order n, we define 
the partition function as the sum of the weights over all elements of C„, 

Zn=Yl ^(^)' (2-2) 

which depends on n, the spectral parameters as well as the functions Xt{-,-)- 
In general, computing Z„ is a hard problem, but something miraculous 
occurs when we set 

Xi[x,y) = X2{x,y) = - 



X3ix,y) = X4{x,y) 



smry 

sin (2 - X + y) (2.3) 



smr/ 

^bix^y) = X6{x,y) = 1, 
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where r/ is a real-valued parameter, known as the crossing parameter. In 
this case, the partition function Z„ = Zn{x; y) can be expressed in terms of 
the so-called Izergin-Korepin determinant [T^, namely as 



sin [xi - yj + ^) sin [xi - yj - ^) 



sm"(" i)(r?) sin(xi-Xj)sin(yi-7/j) 

l<j<j<n 



n 

X det Mix; y), 



where M(x;y) is the n x n square matrix with entries given by 



(2.4) 



sm [xi -yj + j) sm [xi - yj - j) 
The partition function is clearly a trigonometric polynomial, which, by 



(2.4), is seen to be a symmetric function of the Xj's and (independently) of 



the ?/j's. Moreover, it is not difficult to see from (2.2) that it is of degree 
at most n — 1 in each variable. A rather surprising fact about the partition 
function, shown by Stroganov [25j, is that when the crossing parameter r] is 
set at r/ = 27r/3, the partition function becomes completely symmetric in the 
variables {UxjlJUyj}. Stroganov used this insight to find yet another de- 
terminental expression (essentially a Schur function, see \28\ Section 2.5.6]) 
for the partition function, which led to an alternate proof of the Refined 



ASM Theorem and to the formulas (1.3), (1.4) for the "top-bottom" and 
"top-left" doubly refined enumerations. 



From now on, we assume that the weights Xt are given by (2.3) and 
that 7] = 27r/3 (this point in the parameter space of the square ice model 
is sometimes referred to as the combinatorial point). Since the parti- 
tion function is completely symmetric in the x and y variables, following 
Stroganov we use a uniform notation for the spectral parameters, viz. Ui for 
i G {1, . . . , 2n}. To avoid excessive clutter of notation, we reuse Zn = Zn{u) 
to denote the partition function used in this way. We also use the notation 
Z„(ni, . . . , lifc, 0^""'^) for the partition function with the last 2n — k variables 
set to zero. 
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Let 



a(n) = -|,ing+n 

hf \ 2 . \ (2.6) 

b[u) = sm { — — u \ , ^ ^ 



V3 V3 
c(m) = 1, 

and note that the Xt may be expressed in terms of the functions a{-),b{-) 
and c(-) as 

^i{x, y) = X2{x, y) = a{x - y), 

X3{x,y) =X4{x,y) = b{x-y), (2.7) 
Xb{x, y) = Xq{x, y) = l = c{x- y). 

Note that a(0) = 6(0) = c(0) = 1, so in particular when all the spectral 
parameters are set to 0, the weight of every square ice configuration is equal 
to 1, and we get that ^^(0^") = An, the total number of ASMs of order n. 
The following identities are easily proved and will be of use to us below: 

sin(n) = ^{a{u) - b{u)), 

cos(n) = ^(a(ii) + b{u)), 
cos (J^ + u^ = ^(-a(n) + 2b{u)), 
cos (J^-uj = ^(2a(u) - b{u)), 

sin(n — v) = -^{a{u)b{v) — b{u)a{v)), 

sin ( — \- u — v\ = '^{a{u)a{v) + b{u)b{v) — b{u)a{v)). 
V 3 / ^ 

Now, following [23], let 

/„(n) = sin2"-i(n)Z„(n,02"-i), (2.9) 

and note that 

/n(o) = /;(o) = /;'(o) = . . . = = 0, (2.10) 

/(2«-i)(0) = (2n - 1)!Z„(02") = (2n - 1)!A„. (2.11) 



(2.8) 
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The function fn{u) encodes information about the usual "top row" refined 
enumeration of ASMs: by a standard "first row expansion" (dividing square 
ice configurations into classes according to their first row, and noting that in 
the specialization Z„(n, 0^"~^) the contribution to the weight from all other 
rows is constant), we see that it may be written as 



Uu) = sm^^-\u)Y,An,ka''-Hu)b^-\u) 



(2.12) 



k=l 



(see |25l Section 4]). Stroganov also found a fairly simple explicit expansion 
for fn as a trigonometric polynomial: 



Lemma 3. 

fn{u) = K„ ^ 



n— 1 / 4 
n — ^ 



m=Q 



m 



n 



n — m 



sin ((4-3n + 6m)n), (2.13) 



where the constant k„ is given by 



An-, 



(2n-2\ ■ 
\ n-1 J 



(2.14) 



Proof. The formula (2.13) was proved by Stroganov (who used it to obtain a 
new derivation of the formula for A^^k)^ except for the value of the constant 
Kn- To compute Kn, we set u = 7r/3 and evaluate /n(7r/3) in two ways. On 
the one hand, from (2.9) we have that 

n 

Un/S) = sin2«-i (7r/3) ^ A„,,a(7r/3)^-i6(7r/3)"-*^ 



k=l 



V3 



2n-l 



2 /3\" 



^/3 



A. 



n-l 



(since a(7r/3) = l,6(7r/3) = 0); on the other hand, (2.13) gives 



"^^^ / _ 4\ / _ 2 \ 

m=o V / \ / 



n— 1 / 4 
n — o 



m 



R n-i 

<->"-'f-e::;) 



n- 3 
n — m — \ 
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where the last step follows from the Chu-Vandermonde summation identity. 
Comparing the two expressions gives (2.14). □ 

As a side remark, note that one can get a different expression for k„ by 
again equating the two expressions (2.9) and (2.13) and taking the limit as 
ti — 7- (using L'Hopital's rule). This gives 

/n-\ 

Kn = {,-\T-\2n-\)\An 



' n—1 / A 

n — o 



\m=0 



m 



n 



n — m — 1 



(4 — 3n + 6m) 



2n-l 



The fact that this is equal to the right-hand side of (2.14) implies the fol- 
lowing interesting identity. 



Corollary 4. The summation identity 



n 



ni 



m=0 

holds for all n > 1. 



n — m 



(1 — 3n + 6m) 



2n+l 



AY (3n + l)! 
3 



(2.15) 



An independent proof of (2.15) was found by Amdeberhan [T]. 

Our next result will give a hint that the remarkable function fn{u) also 
holds the key to higher order refined enumerations of ASMs, which are 
related to the A;-variate specialization Z„(ui, . . . , n^, 0^"'"'^). 

Theorem 5 (Partition function formulas). The partition function with all 
but k variables set to zero can be written as 



Zn{ui, ... ,Uk,0 



2n-k\ 



Cn,k 



dui ^ 



(2.16) 



JJsin^" ''{ui) Yl sm{ui-Uj 

i=l ^<i<j<k 



where the constant C,n.k is given by 

Cn,k 



^n"'ni=i'(2n-jr-^- 
In particular, taking k = 2n we get the expression 

1 



(2.17) 



Zn{ui, ■ ■ ■ ,U2n) 



det ' ^ ' 



dui 



I n 

l<i<i<2n 



(2.18) 

for the full partition function, where H{n) = 13^=1 ^'^ the hyperfactorial. 



14 



We note in passing that several different representations for the partition 
function have been found, see, e.g., fT^ [71 1^ IT7 1 1^. 
The proof wiU require the following lemma. 

Lemma 6. Let ri, . . . , r^v be distinct real numbers, and let < k < N . The 
solutions to the homogeneous system of linear equations 

N 

Yl ^h.-oJt = (1 < J2, . . . , Jfc < iV, < m < iV - A: - 1), (2.19) 
ii=i 

^'^(j-O.-.^Ofc) = sgn(fT)zji,...j, (cre5fc), (2.20) 

in the family of unknowns (-Zji,...jfc)i<ji,...,jfc<Af) ^'"e determined up to a global 
multiplicative constant. 



Proof. The anti-symmetry condition (2.20) implies that any variable Zj^ 



where j'l, . . . ,jk are not all distinct is 0. It will therefore be enough to show 
that for any two index vectors (ji, ■ ■ ■ , jk), U'l^ ■ ■ ■ J'k) where ji, . . . ,jk are 
distinct and j'l, . ■ ■ ,j'f^ are distinct, the ratio Zji^^ j'^ is determined 

uniquely by the equations. Furthermore, if we prove this claim in the case 
that the vectors (ji, . . . ,jk), (ji, • • • ,jk) differ in only one position then the 
more general claim follows by successively replacing one index at a time 
(using anti-symmetry to permute the order of the indices). Again because of 
ant i- symmetry, it suffices to show this when the vectors differ only in the first 
position. That is, we wish to determine the ratio -2ji,...,jfc/2^jjj2,.--jfc '^here 



jij j[: 32, ■ ■ ■ , jk are distinct. Note that the equations (2.19) are equivalent 
to the assumption that for any polynomial p{x) of degree < — A; — 1 we 
have the relation 

N 

J2^»d2,...,hPira) = 0. (2.21) 

a=l 

We apply this for the polynomial 



1</3<]V 



In this case, all the summands in (2.21) vanish except when a = ji,ji, so 
we get that 

P(^ii)^iij2,...jfe +Pirj[)zj'^j^,...,j^ =0 
(where we know that p(rj-^),p(rji^) ^ 0), which determines the ratio of the 
two variables Zjx,j2,---dk/ ^j[j2,-,jk-' claimed. □ 
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Proof of Theorem^ Define the function gn^k as 

gn,k{ui, ...,Uk) = Zniui, . . . ,14^,0^"" 

X JJsin^"'~''(ui) Y\ sm{ui-Uj). ^ ^ 

i=l ^<i<j<k 

The proof is based on the observation that Qn^k satisfies the following list of 
properties: 

1. Qn^k is a trigonometric polynomial in the Uj's of degree at most 3n — 2 
in each variable. 

2. gn,k{---,Ui + 7r,...) = {-lYgn,k{- ■ ■ ,Ui, . . .). 

3. gn,k{- . . ,Mi, . . . ) +gn,k{- . . ,ni + ^, . . . ) + gn,k{- . . + ^, . . . ) = 0- 

4. gn ]^ is divisible by sin^"~^(ui) for each i. 

5- 9n,k{- • • 1 '^ij ■ ■ ■ , Uj , . . .) = gn,k{- • ■ ) ^ji • • • j ^i, . . . ). 

Property 1 follows from the fact that the partition function Z„ is a 
trigonometric polynomial of degree < n — 1 in each variable. Property 2 
follows from the relation 

Z4...,u^ + 7T,...) = (-l)"-iZ„(. ..,Ui,...) (2.23) 

satisfied by Z„, which is a consequence of the fact that a{u + vr) = —a{u), 
b{u + vr) = —b{u), together with the observation that each row and column 
of a square ice configuration contain an odd number of vertices with weight 
c (translate this to the language of ASMs to see why). Property 3 is a 
consequence of the symmetry of Z„ in the variables ui, . . . , U2n together 
with the following functional equation proved by Stroganov [25]: 

Fniui, ...) + F„ (m + f ,...)+ F„ (m + f ,...)= 0, (2.24) 

where Fn{ui, . . . , U2n) is defined by 

2n 

Fn{ui,...,U2n) = (^1 , • • • , ti2n) JJ sin(Mi - Uj). 

3=2 

Properties 4 and 5 are immediate from the definition of gn,k- 

We will now show that any function which satisfies these five properties 
is determined uniquely up to a normalization constant depending only on n 
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and k. The idea is a generalization of Stroganov's approach to the proof of 
the Refined ASM theorem [25J. Starting from the first property, we express 
Qn k in canonical form, 



3n-2 



3n-2 



9n,k\Ul, 



ai=— 3n+2 afc=— 3n+2 



Cai,...,ak 6Xp 



(2.25) 



where the Fourier coefficients Cai,...,afc are to be determined. 



Using the second property, we see that Cai,...,afc = whenever some aj 



does not have the same parity as n. Using the third property, Ca^, 
whenever some is divisible by 3. So we may rewrite (2.25) as 



where 

/„, = {2 — 3n < a < 3n — 2 : a = n (mod 2), a = 1,2 (mod 3)}. 
Note that |/n| = 2n. Now, the fourth property implies that 

Qrr 



duf 



gn^kiui, ■ ■ ■ ,Uk) = 0, {0 < m <2n - k - 1), 



ui=0 



which translates to the equations 

Cai,...,akaT = (02, ..■,ak £ In, 0<m<2n-k-l) 



(2.26) 



on the coefficients. Similarly, the anti-symmetry property 5 implies the 
corresponding anti-symmetry condition 

Caia^),...,a{ak) = Sgn(cT)Cai,...,afe (2.27) 

on the coefficients, which holds for any permutation a of the elements of 
/„. But we now see that (2.26) and (2.27) are identical to the system of 
equations (2.19), (2.20) from Lemma [6| with the substitution N = 2n and 
with ri, . . . ,rjv denoting the elements of I„. The conclusion of Lemma [6] 
guarantees that the coefficients Cai,...,afe5 and hence the function gn,k-, are 
determined uniquely up to a multiplicative constant, which was our claim 
up to this point. 
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Next, we claim that 



hn,k{ui, ■ ..,Uk) 




-1 



(2.28) 



satisfies the same hst of properties discussed above. First of ah, fn{u) is a 
trigonometric polynomial of degree 3n — 2 in u and taking its derivatives 
cannot increase the degree, so /i„ satisfies the first property. The second 
and third properties follow immediately from the equations 

/„(n + vr) = (-l)"/„(tx), 

= fniu) + + y) + fn (u + y 



satisfied by fn{u) (the first equation is again a consequence of ( 2.23| ), and 
the second equation is the special case U2 = ■ ■ ■ = U2n = of (2.24)). 



The fourth property holds because the smallest factor of sin(u) in a 
derivative of /„ is sin^"~'^(ii) in the {k — l)th derivative. The fifth and 
last property is immediate. Thus, we have shown that gn{ui, . . . ,Uk) = 
Cn,khn{ui, ■ ■ ■ ,Uk) where Cn,k is a proportionality constant. 

It remains to prove the formula (2.17) for Cnk- We shall do so by an 



induction on k for each fixed n. We first note that Cn,i = 1 by (2.9). We 
now evaluate Cn k by taking the limit of all Uj's to zero in (2.16), 



Cn,k 
An. 



n 



sm 



2n~k 



(ui) sin(nj 



lim 

Ul,...,Uk- 



l<i<j<k 




(2.29) 



We compute the limit as an iterative limit, by taking the variables succes- 
sively to 0; in fact, because of the inductive nature of the computation we 
only need to let — )• since a recursive structure is revealed. To do this, 
we use L'Hopital's rule by taking the (2n — k)th partial derivative of the nu- 
merator and denominator with respect to Uk and setting each to 0. Consider 
first the relevant terms in the numerator. We have to evaluate 



d' 



,2n~k 



k-1 



sin^" ''{uk)Ylsm{ui 



Uk) 



(2.30) 



This multiple derivative will lead to a sum involving multinomial coefficients 
over all possible combinations of derivatives of the k terms inside in general. 
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However, since we are going to set Uk = at the very end, the only term 
that contributes is the one where ah 2n — k derivatives act on sin^"~'^(ufc). 
We thus obtain, after easy manipulations. 



(2n - A;)! J| sin(tii). 



(2.31) 



i=l 



We now consider the denominator. Since Uk is only present in the last row of 
the matrix, we can perform a cofactor expansion of the determinant using 



the last row. By (2.10), the only term that contributes after taking the 



(2n — k)th derivative is the cofactor of the entry in position {k,k), which 
leads to 



2n-l 



Uk=0 



i2n-l)\Andel I flllp^^ 
«'i=i V dui 



(using (2.11)). To summarize, we have shown that 

fc-i 



Cn,k ^ {2n-ky. 

An {2n — ly.An ui,...,Uk-i^O 



lim 



(2n-fc)! Cn,fc-i 
(2n-l)!^n An 



n 

i=l 



sm{ui-Uj) 

i<t<j<fe-i 



(2n-l)!A„det ( '^''{"1"'^ 

i,j=l \ dul 



(2.32) 



It is easy to check that the expression on the right-hand side of (2.17) satisfies 
this. □ 



When one is interested in refined enumeration formulas, one would like 
to compute Zn{ui, . . . , Uk, 0^"~^) for large n and fixed k. In this case (2.16) 
is very useful because the order of the determinant is fixed. Note, however, 
that this formula would not be very practical to use if the computation of 
the derivatives of /„ were not efficient enough. As we shall see in the next 
section, for the purpose of deriving enumeration formulas we are interested 
in a representation of the derivatives /n™^(^i) not as a standard trigono- 



metric polynomial (this would be immediate from (2.13)), but instead we 
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need a specific expansion involving products of powers of a{u) and b{u) fac- 
tors, similarly to (2.12), with an additional factor consisting of a power of 
sin(M). We now show that there is an efficient recursive algorithm for the 
computation of such a representation. 

Proposition 7. For < m < 2n—l, the mth derivative of fn can he written 
as 



firHu) = sin2"-i-™(n) ^ c„,„,fca^-^(n)6"-'=+-(n), 

k=l 

where c„^o,fc = ^n,k) f^iT-d for m > we have the recurrence 

Cn,m,k =5 ( ~ 2A;C„,rrt-l,A;+l + {n + Ak - 2m - 2)Cn,m-l,k 

+ (5n - 4/c + 2m + 2)cn,m-i,k-i - 2{n - k + m + l)cn,m-i,k-2^ ■ 

(2.34) 

Proof. Since the argument of a and b will always be u, we will omit it in 
this proof. Starting with 



n+m 



(2.33) 



/(™-i)(n) = sin^"— (n) ^c„,^_i,fea'=-i6"+'"-i^ 



(2.35) 



and differentiating, we get 



/(™)(^x) =(2n - m) sin2"-i— (n) cos(u) ^ 



.1 fca'=-^6'"+™-i-'= 



+ sin2'^ "'(n)^c„,m_i,fc 



(fe - l)a'=-2^"+— i-'^ A cos + u 
V 3 ^3 



(n + m - 1 - k)a^'^h 



fc— li n+m— 2— fc 



COS u 

^/3 V3 



vr 



(2.36) 



Now we take out the common factor of sin ~"^{u) and use the identities 
(2.8) to replace all trigonometric factors by sums of products of a's and 6's. 
After simplifying and collecting terms, we get 



f!r\u) =^sin2"-i-™(t,)^c„,^_i,, 

k 

-,A;— lin+m— fc 



2{k - l)a'''^n-k+m+l 



+ {n-2m-2 + ik)a^'^b'' 
+ (5n + 2m - 2 - 4A:)a^6"+'"-i-'= 
- 2(n + m - 1 - fc)a'=+^6"+™-2-'= 



(2.37) 
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which is another way of encoding (2.34). Note that if we assume inductively 
that the range of values of k for which the coefficient Cn,m-i,fc is nonzero is 
l<fc<n + m — 1, then the range of fe's over which we should consider 
Cn,m,k isO<A;<n + m + l. However, the coefficient 2k appearing in front 
of the factor Cn,m-i,k+i in (2.34) ensures that c„_m,o = 0; and similarly the 
coefficient 2(n — k + m + 1) of c„^m-i,fc-2 ensures that Cn,m,n+m+i = 0, so 
in fact we only get nonzero coefficients Cn,m,k when 1 < k < n + m. This 
explains the range of summation in (2.33). □ 



The proof of Proposition [7] suggests a new way of looking at the partition 
function Zn- The usual way of considering it as a trigonometric polynomial 
has been replaced by a polynomial in a and b. This is in some sense simpler 
because the total degree of the polynomial is constant. Underlying this is 
a change-of-basis transformation in the space of trigonometric polynomials 
which has nice properties. This point of view deserves to be examined in 
more detail. In particular, one interesting open problem would be to find 
closed (non-recursive) formulas for the coefficients Cn,m,k- 



Note also that the recurrence (2.34) is not quite as arbitrary as it might 
appear at first sight. The coefficients of Cn,m-i,k on the right-hand side obey 
a symmetry relation under the interchange k-^n + m + l — k, namely 

-2k o -2{n -k + m + 1), {n + Ak - 2m - 2) {5n - Ak + 2m + 2). 

The significance of this observation is unclear to us at the moment. 



3 Boundary behavior: triple and quadruple re- 
finements 

We will now prove Theorems [T] and [2] The idea is to equate the two formulas 



for the specialized partition function, one given by the definition in (2.2), 



and the other by the determinant in (^2.16) 



We now set the notation for use in the proofs. As mentioned before, we 
will consider f^\u) to be a polynomial in a{u) and h{u). This polynomial 
turns out to be homogenous of degree 3n — 2. Therefore, it makes sense to 
consider it as a polynomial in the single variable t = a{u)/b{u). For conve- 
nience, we will use the same notation f^\t) for the new function (where 
the derivative is still with respect to the variable u) . From Proposition [7] 
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and from (2.8), it follows easily that 




2n—l—m 



uin-2 



{U) (t - 1) 



2n— 1— m 



n+m 
k=l 



fc-1 



(3.1) 



We will have at most four variables in the determinant we want to evaluate. 
In continuation with Theorem [5} we use Ui for the variables. Throughout 
the computation, let Oj (resp. bi) be a shorthand for a{ui) (resp. b{ui)), and 
denote ti = ai/hi. 



3.1 Triple refinement 

Proof of Theorem [7| We will use Theorem [5] to calculate the number of 
ASMs where the 1 in the first row is in the ith column, the 1 in the first 
column is in the jth row, and the 1 in the last row is in the A:th column. To 



relate the determinant in (2.16) to the generating function of the numbers 



j4^^^(i, j, /c), we need to calculate a determinant of a matrix of size 3 of 
derivatives f^\ud), corresponding to m = 0,1,2, expressing each of the 



derivatives in the form (3.1). Let us compute these derivatives; we shall see 



that the result will involve the functions Q!„(t), /3„(i), 7n(t) from (1.5). For 
m = we have 



V3 



2n-l 



k=l 



2n-l 



2n-l 



Something interesting happens while computing fn\t). Applying the re^ 



currence (2.34), we get 

2n-2 




L3n-2 



{u)it-l 



>2n-2 



2kA^,k+i + (n + 4/c - 4:)An,k 



(3.2) 



k=l 



+ (5n - 4fc + 4)^„,fc_i - 2{n-k + 2) A„_fc-2 1 1 



Lk-l 
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Curiously, the coefficient of ^ in tliis expression simplifies to a constant 
(depending on n but not on k) times A^-i^k-i^ because of the identity 



(n-2)!(3n-2)! 
(2n-3)!(2n-2)! 



An-i,k-i = - 2kAn,k+i + {n + 4k- 4) A„,fc 

+ (5n -4k + 4)An,k~i -2{n-k + 2)An,k-2 



(which can be easily verified from (1.2) by elementary algebra). Thus we 
get that 



n-l 



2)!(3n 2)! 3„_2 2n-2^ (3 3) 

4/ 2(2n-3)!(2n-2)! ^ ^ ^ 



Unfortunately, this nice phenomenon does not persist for higher deriva- 
tives. For the second derivative, using (2.34) again, after a short computa- 
tion we get 



V3' 

n+2 



2n-3 



(n-2)!(3n-2)! 3„_2^^^^^_^^2n-3 



4(2n-3)!(2n-2)!' 
^(-2{n-k + 3)An-i,k~3 + (5n - Ak + 6)An-i,k-2 



k=l 



(3.4) 



+ {n + 4k- 6)A„_i,fc_i - 2M„„i,fc t 



,k-l 



2n-3 



(n-2)!(3n-2)! 3„_,^^^^^_^^,„_3^^^^^^ 



4(2n-3)!(2n-2)! 



(The expression for 7n(i) does not seem to simplify further.) 

Let us now look at the denominator on the right-hand side of (2.16). In 
the case k = 3, it can be written (making use of (2.8) again) as 



n 

i=l 



sin2"-3(ui 



sm{ui-Uj) 
i<j<i<3 

3n-3 3 



i=l 
3n-3 3 



Jj52n-3(^._l)2n-3 "Q b,b,{U - t,] 

1<«<J<3 



i<«<i<3 



i=l 
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Combining this with the above computations, we find that in the case k = 3, 



(2.16) transforms into 

^n(ni,n2,n3,02"-3) 



X det 



ih-lfanih) {t2-lfan{t2) - if Units) 

ih - l)/3n(tl) (t2 - l)/3n(t2) (^3 " l)/3n(t3) 
7n{tl) 7n{t2) ln{h) 



(3.5) 



where pn = |Cn,3 
the value of Cn3 from (|2.17l)). 



(n-2)!(3n-2)! 
(2n-3)!(2w-2)! 



(which becomes (1.7) upon substituting 



Now we move on to the left hand side of ( 2.16 ) with = 3. Our goal is to 
use the original definition of Zn{ui,U2, u^, 0^"""^) to represent it as a factor 
of 6162^3 times a generating function in ti,t2,t3 expressed in terms of the 
coefficients A]^^^. Equating this expression to (3.5) will give the theorem. 



First, note that since Z„ is a symmetric function, we are free to choose 
which of the spectral parameters xi, . . . , Xn, yi, ■ ■ ■ ,yn each of the variables 
ui,U2,U3 is chosen to stand for. Since we are interested in relating Zn to the 
coefficients the natural choice is to set xi = ui, yi = U2 and Xn = U3; 

all other spectral parameters are of course set to 0. For this choice, the 



weight w{C) of each square ice configuration in the sum (2.2) is influenced 
only by the behavior of the configuration in the top and bottom rows and 
the leftmost column. In particular, we can divide the configurations into 
classes according to the position of the Is in the top and bottom rows and 
leftmost column in the corresponding alternating sign matrix. Within each 
class, all configurations have the same weight. 

Let C G Cn he a square ice configuration, and denote by k the po- 
sitions of the unique 1 in the top row, leftmost column and bottom row, 
respectively, of the corresponding ASM. There are three possible cases, each 
giving rise to a separate computation of the weight w{C). 

1) i = j = 1, 2 < k < n. The number of configurations of this type is 
^n^^(l; 1; ^) = ^n-i,fc-i> and the weight associated to such a configura- 
tion is given by 



n— 1 „n— 2 n—kik—2 



a 



an 



^'b{u3 - U2) 



n— 1 2 n~kik~2 



a 



(3.6) 



(&lft2&3) 



n-1. 



1 — ^3 + ^2*3)^2 ^^3 



These weights arise as follows: in the first row, the first vertex has weight 
c(ui — U2) = 1, and the remaining n — 1 vertices have weight h{ui) = hi. 
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In the first column, the vertices in positions 2 through n — 1 have weight 
b{—U2) = a{u2) = 02, and the last vertex has weight b{u3 — U2), which 
can be rewritten as a2a3 + 62^3 — ^2^3 (using (2.8)). In the last row, the 
vertices in positions 2 through k — 1 have weight b{u3) = 63, the vertex 
in position k has weight c(m3) = 1, and the vertices in positions k + 1 
through n have weight 0(^3) = 03. 

2) 2 < i < n, j = n, k = 1. The number of configurations is and 
the weight is seen (using similar reasoning to that explained above) to 
be 

w{C) = a\-^'l-'a]^-^b^-'^a{ui - U2) 

= a\-%'l-'a'^-^bl^-\aia2 + 6162 - ^laa) (3.7) 
= {bib2bsr-\l-t2 + ht2)t\-H^-\ 

3) 2 < i,k < n, 2 < j < n — 1. In this case the number of configura- 
tions is A'^^^{i, j, k). A short computation gives the weight of such a 
configuration as 

w{C) = a^^h^T'^f^-^h^^-^al-Hl-^aiui - 1^2)6(^3 - U2) 
X (aia2 + 61&2 - bia2){a2a3 + 62^3 - ^203) 

= (6l62&3)"-'(l -t2+ tlt2){l -h + t2h)t\-X~'-H^-\ 

Now all that is left to do is to combine the three cases to arrive at the 
partition function. We have 

n 

— — — ^Z„(ni, n2, n3, O^"-^) = {1 - + ^2*3) V A^^i^k-it^^t^''" 
(010203) ^ 

n 



1=2 

n n—1 n 

{i-t2+ ht2){i - 13 + t2t3) E E E m-Hf'tr". 

1=2 j=2 k=2 



Comparing this to the result (3.5) gives us just the equality of generat- 
ing functions claimed in (1.6) (except with the variables x,y,z replaced by 
^1,^2,^3), and finishes the proof. □ 
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3.2 Quadruple refinement 



Proof of Theorem [1[ The strategy is the same as that of Theorem [TJ Recall 
that 

An^^^{i,j,kJ) = #\^{mij)lj^i G ASM„ | mi,i = mj,i =mn,fc = m£,„ = l}. 

(3.8) 

As before, we will use Proposition [7 to calculate fn^ . Since the sum- 

(2) yl 

mand in the expansion of /A (t) in ( |3.4 ) is quite complicated, we leave it 
to the reader to verify that the recursive formula indeed yield a function 
proportional to 6n{t) ( |1.5[ ), 

2n-4 



n+3 



fc-1 



fc=l 



n-2 



(^_2)!(3n-2)! ^3,. 
8(2n-3)!(2n-2)! 



(3.9) 



^(n)(t-l)"5„(t) 



The denominator on the right-hand side of (2.16) in the case k = A can 
be written as 



n 

i=l 



sin2"-4(^i,- 



n 

l<i<j<4 



4n-5 4 



n^f"'(ti-i)" n 

i=l l<«<j<4 
4n-5 4 

\{hi--\u-if--^ n (*^-*^)- 

i=l l<j<j<4 



/n\ 

Combining the calculation of /„ (t) in ( 3.9 ) with those in ( 3.3 ) and ( 3.4 ) , 
we find that in the case /c = 4, (2.16) transforms into 

Z„(U1,U2,U3,U4,02"-^) = C7„(6i62&3b4)""^A(ti,t2,i3,i4)"^ 

/ (ti-l)3a„(ti) (t2-l)='a„(t2) {t3-lfan(k'.) iU-iy^aniU) \ 

(tl-l)2/3„{ti) {t2-l)2/3„(t2) (t3-l)^/3n(t3) {U-l)"^ P„{U) 

(tl-l)7n(tl) (i2-l)7n(i2) (t3-l)7n(t3) (t4-l)7n{t4) 

V <5n(tl) <5n(t2) <5n(t3) 5n(t4) / 

(3.10) 



X det 



where (T„ 

of (fTTol). 



1 Cn,' 



(n-2)!(3n-2)! 
4 I (2n-3)!{2n-2)! i ' 



which is the same as the right-hand side 
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Now we compute the partition function Z{ui,U2,us,U4,0^"'~^) from the 
definitions. This time we will set the spectral parameters as Xi = U\,y\ = 
U2,Xn = and ?/« = ^4- Let C € C„ be a square ice configuration, and we 
will use the variables i, j, k, t to denote the positions of the c-type vertices in 
the first row, leftmost column, last row and rightmost column, respectively. 
These correspond to the spectral parameters u\,ui-,U2,-,U4. Remember that 
i and k are column parameters and j and I are row parameters. Now there 
are seven possible cases depending on the values of i, j, fc, I. 

1. i = J = 1, /c = f = n. The number of such configurations is clearly 
An^2- The weight is given as follows: both the (1,1) and the (n, n) 
vertices get a weight of c{ui — U2) = 1 and c(u3 — U4) = 1 respec- 
tively. All other boundary vertices are of 6-type. The other two 
corner vertices at (l,n) and (n, 1) get a weight of b{u3 — U2) and 
b{ui — U4) respectively. Therefore we have n — 2 factors of each of 
b{ui) = W2) = a2,b{u3) = 63 and 6(— U4) = 04. Thus the total 
contribution to Zn from such configurations is 

^n-2(&l 02^304)"^"^ b{u3 - U2) b{ui - U4), 

= At-2(^ia2^3a4)"~^(a2a3 + 62^3 - 6203) (01^4 + ^1^4 ~ Mi), 

= An-2{blb2b3b4r-^ (1 - is + t2t3)(l - + ^4^1) (^2*4)''"^ • 

(3.11) 

2. i = j = n,k = I = I. The reasoning in this case is essentially identical 
to that of the case above. This time all boundary vertices except the 
ones at (l,n) and (n, 1) are of a-type. We eventually obtain a total 
contribution of 

An-2ibib2b3b4T'^ (1 - i2 + iii2)(l -U + hU) (iiis)""' • (3.12) 

3. i = j = 1,1 < k,£ < n. The next four cases follow the same kind of 

reasoning as the one here. We describe this in detail and leave it to 
the reader to verify the next three cases by analogous reasoning. The 
number of configurations of this type is 

Al^^^il,l,k,i) = A^^,{k-l,i-l), 

which is the number of ASMs of size n—1 with a 1 in the bottom row 
at column A; — 1 and a 1 in the rightmost column at row i — 1. The 
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weight associated to such a configuration is given by 

=6r2ar'ar'"''^r'«l~'^r'''(«2a3 + ^263 - ^2^3) 
X {aia^ + 6164 — 640i)(a3a4 + 63^4 — 6304) 

= {blb2b3b^r-Hr^t^-'-%-\l -ts+ t2t3) 
X {l-t4 + ht4)il-ti+Uti). 

(3.13) 

These weights arise as follows: in the first row, the first vertex has 
weight c{ui — U2) = 1, the next remaining n — 1 vertices have weight 
b(ui) = bi and the last vertex has weight ^4). In the first column, 
the vertices in positions 2 through n — 1 have weight 6(— ^^2) = cl{u2) = 
02, and the last vertex has weight b{u^ — U2), which can be rewritten as 



0203+62^3 — ^203 (using (2.8)). In the last row, the vertices in positions 
2 through k — 1 have weight 6(ti3) = 63, the vertex in position k has 
weight c{u^) = 1, the vertices in positions k + 1 through n — 1 have 
weight 0(7x3) = 03, and the last vertex has weight a(n3 — 114). In 
the last column, the vertices in positions 2 through i — 1 have weight 
6(— U4) = 04, the vertex in position i has weight 0(^4) = 1 and the 
vertices in positions £ + 1 through n — 1 have weight a(— U4) = 64. 

We now note that A^^^{k -1,^-1) = Al^^{n + 1 - k,n + 1 - £) hy 



rotational symmetry. Using the generating function defined in (1.8), 
the sum over k,£ after elementary manipulations is given by 

(M2M4)"-' tr' t2-H^-t3+t2t3) 

il-U + -h+ Ml)^r^-l (^t3, ^ 

4. j = n,k = 1,1 < i,£ <n. There are A^'^^{i-1,£) = A^^^{n + l-i,£) 
such configurations and the weight of such a configuration is given by 



W{C) ={bib2b3b,r-H\-^ 

X (1 - t2 + ht2){l -t4 + t3U){l -ti+ Uh) 

The sum over i,£ is now given by 

(M2M4)"'' ^r' t'i-Hi-t2+tit2) 

il-U + t3t4)(l -ti+ Uti)Alh^ (u, 



(3.14) 



ti 
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(3.15) 



5. k = £ = n,l < i,j < n. There are A'^}^-^{i,j) such configurations and 
the weight of such a configuration is given by 

w{c) ={bib2b3hr-H\-'' tr'~' 

X (1 - i2 + ht2){l -t3+ t2t3){l -ti+ t4tl) 

The sum over i,j is given by 

(61626364)""' *r^(i - *3 + ^2*3) 

(1 - i2 + tii2)(l -ti + Uti)Al}l^ (h, ^ 



t2 

6. i = n,e=l,l <j,k <n. There are A]^^;i(fc,j-1) = A^}:^{k,n+l-j) 
such configurations and the weight of such a configuration is given by 

w{c) ={bib2bshr-Hr' tr'-' tr'-' 

X{l-t2+ ili2)(l - ^3 + t2t3){l -t4 + t3U)tl-'^. 

The sum over j, k is given by 

(61626364)"-'*?-' tr'(l-t3 + t2t3) 

(1 - ^4 + t3U){l -t2 + tlt2)Albi (^t2, . 

7. 1 < k,£ < n. The number of such matrices A'^^^^{i,j, k,£) is the 
quantity we are ultimately interested in calculating. We describe in 

detail the weight of vertices in the first row. The weights of vertices 
in the first column, last row and last column are obtained by similar 
arguments. In the first row, the first vertex has weight a{ui — U2), the 
next i — 2 vertices have weight a{ui) = oi, the ith vertex has weight 

c{ui) = 1, the next n — i — 1 vertices have weight b{ui) = 61, and the 
last vertex has weight b{ui — U4). The weight comes out to 

w{C) = a\-^ Vl-'-^ al-^-^ 6^' ^-^^ 6,^-^ V^-x-i 

X o(«i — U2)b{u3 — ^2)6(1*1 — U4^)a{u3 — Ui) 

X (0203 + 6263 - 6203) (aia4 + 6164 - b^ai) 
X (0304 + 6364 - 6304) (aia2 + 6162 - 6102) 
=(61626364)"-' q-i-j t^-i-kte-2 

X{l-t2+ tit2){l -t3 + i2*3)(l - ^4 + t3t4){l - h + Utl)- 

(3.17) 
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As in the triply-refined case, comparing the sum of the contributions from 
these seven cases to the result (3.10) gives us the claimed identity (1.9). □ 
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